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Abstract 



Oh. We consider the fixed-dimension perturbative expansion. We discuss the nonanalytic- 

c"| ! ity of the renormalization-group functions at the fixed point and its consequences for the 

numerical determination of critical quantities. 

In the last thirty years there has been a significant progress in the understanding of critical 
phenomena. It has been realized that the behavior in the neighborhood of a critical phase 
transition, i.e. a transition characterized by long-range correlations, is determined by very few 
properties, the space dimensionality, the range of the interactions, the number of components of 
the order parameter, and the symmetry of the Hamiltonian. This means that physically different 
systems may have the same critical behavior. For instance, a simple fluid at the liquid-vapor 
transition and a uniaxial magnet at the Curie point behave identically: Critical exponents, 
dimensionless amplitude ratios, scaling functions are numerically equal. This phenomenon, 
which is referred to as universality, has been understood within the Wilson's renormalization- 
group (RG) approach. The conceptual setting is thus quite well established, and the theory 
of critical phenomena has thus reached the maturity of well-verified theories like, for instance, 
QED or the standard model of weak interactions. Nonetheless, it is important to improve 
experiments and theoretical calculations in order to understand the limits of validity of these 
theories. In order to test QED and the standard model, several experiments have provided 
accurate estimates that can be directly compared with the theoretical predictions. The most 
classical ones are the experiments on the ^-factor of electrons and muons, and on the Lamb shift 
in hydrogen. In the theory of critical phenomena, the superfluid transition in 4 He plays a very 
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special role, since it is essentially the only case in which one can determine a critical exponent 
with an accuracy of 10~ 4 . This is due to a combination of particularly favorable conditions: the 
singularity in the compressibility of the fluid is particularly weak; it is possible to prepare very 
pure samples; experiments may be performed in a microgravity environment (on the Space 
Shuttle, for instance), thereby reducing the gravity-induced broadening of the transition. A 
recent experiment |], || obtained an extremely accurate estimate^] of a, a = —0.01056(38). 
This result should be compared with the most precise theoretical estimates: The analysis of 
high-temperature (HT) expansions gives|, [§ -0.0146(8), -0.0150(17); Monte Carlo (MC) 
simulations give|| [f| —0.0148(15), —0.0169(33); The d = 3 perturbative expansion gives || |7j 
—0.0112(21), —0.011(4). There is a clear discrepancy between the most accurate theoretical 
estimates and the experimental result. However, in order to understand if the difference is truly 
significant, we must ask the question: 'Are the quoted errors reliable?' Our experience, looking 
backward in time, is that there is a natural tendency to be overconfindent in own's results, 
and thus to systematically underestimate the errors: As Hagen Kleinert || put it, each one 
has a tendency to apply the "principle of maximal optimism." Clearly, further theoretical and 
experimental investigation is needed to settle the problem. 

In order to set correct error bars, it is necessary, although clearly not sufficient, to have a 
good understanding of the possible sources of systematic error. In MC and HT works, most 
of the systematic error is due to the nonanalytic corrections to scaling. Indeed, in iV-vector 
systems, there are corrections t A (t is the reduced temperature) to the leading scaling behavior, 
with A pa 0.5-0.6 in the physically relevant cases < iV < 4. In the analysis of the MC data 
and of HT series, these slowly decaying corrections require careful extrapolations, in the absence 
of which precise but incorrect results are obtained. To give an example, we report here some old 
(but not too old!) results for the four-point renormalized coupling g* in the three-dimensional 
Ising model (see the discussion in Sec. 5 of Ref. || and Fig. 1 reported there): 



MC, no nonanalytic corrections: |J 
MC, with nonanalytic corrections: [Q 
HT, no nonanalytic corrections: [TO] 
HT, with nonanalytic corrections: [|TT 



g* = 25.0(5); 
g* = 23.7(2); 
g* = 24.5(2); 

g* = 23.69(10), 23.55(15). 



For comparison, perturbative field theory gives g* = 23.64(7), while a recent analysis of 
improved HT expansions gives |12| g* = 23.49(4). Clearly, neglecting the corrections to scaling 
introduces a large systematic error! And, even worse, there is no way to evaluate it, unless one 
assumes that nonanalytic corrections are really there. 



A solution to these problems is represented by the so-called improved models, [13, 14, 15 



IB] , |5], [T2|, |], [3| which are such that the leading scaling correction (approximately) vanishes. The 
systematic errors are now sensibly reduced and one obtains more reliable estimates. 

MC and HT analyses, although different in practice, are very similar in spirit, and indeed 
they are affected by the same type of systematic errors. In order to assess the reliability of the 



lr The original result reported in Ref. El was incorrect. The new estimate is reported in Ref. The error 
is a private communication quoted in Ref. Hi. 
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results, it is thus important to have a completely different approach to compare with. Field 
theory provides it and indeed independent estimates can be obtained using a variety of differ- 
ent methods: the e-expansion pioneered by Wilson and Fisher, the fixed-dimension expansion 
proposed by Parisi, the perturbative expansion in the minimal-subtraction scheme without e- 
expansion proposed by Dohm, and the so-called exact RG, which essentially consists (there are 
many different versions, see, e.g., Ref. [17|) in approximately solving nonperturbatively the RG 
equations. 

In this contribution we will focus on the fixed-dimension expansion method, which is the 
one which has provided the most precise estimates, and, together with the e-expansion, has 
been the most used. We consider the standard 4 theory with iV-vector fields and discuss the 
role of the singularities of the RG functions at the critical point in the numerical determination 
of critical quantities. 



1 Singularities of the RG functions 



An important controversial issue |y|, [19], ^(| [H], |22| in the field-theoretic (FT) approach in 
fixed dimension is the presence of nonanalyticities at the fixed point g*, which is defined as the 



zero of the /3-function. The question was clarified long ago by Nickel |19j who gave a simple 
argument to show that nonanalytic terms should in principle be present in the /5-function. The 
same argument applies also to other series, like those defining the critical exponents: Any RG 
function is expected to be singular at the fixed point. 

To understand the problem, let us consider the four-point renormalized coupling g as a 
function of the temperature T. For T — > T c we can write down an expansion of the form 

l + ait + a 2 t 2 + ... + M A + ht 2A + ...+ 

Cl t A+l + ... + d^ 2 + ... + e x V + ..], (1) 



9 = 9 



where A, A 2 , . . . are subleading exponents and t is the reduced temperature t = (T — T c )/T c . 
The corrections proportional to V are due to the presence of an analytic background in the 
free energy. We expect on general grounds that a\ — a 2 — a s — . . . — 0. Indeed, these analytic 
corrections arise from the nonlinearity of the scaling fields and their effect can be eliminated 



in the Green's functions by an appropriate change of variables. |23| For dimensionless RG- 
invariant quantities such as g, the leading term is universal and therefore independent of the 
scaling fields, so that no analytic term can be generated. Analytic correction factors to the 
singular correction terms are generally present, and therefore the constants q in Eq. ([!]) are 
expected to be nonzero. 

Starting from Eq. ([!]) it is easy to compute the /3-function. Since the mass gap m scales 
analogously, for A < 7 (this condition is usually, but not always, satisfied^), we obtain the 
following expansion: 

^9) = m^- = ai Ag + a 2 (Ag) 2 + . . . + ^(Ag)* + /3 2 (Ag)* + ...+ 

7 i(As) 1+ * + • • • + d^Ag)^ + ... + Ci(Ag)i + . . . , ( 2 ) 



2 In some models, for instance in the 2D Ising model, A > 7. In this case, Eq. (||) is still correct [^i) if 7 
and A are interchanged. Also u\ — —j/v in this case. 



3 




Figure 1: Distribution of the results for the resummations of g* for p = 1/10, c = —Zi/5. 



where Ag = g* — g. It is easy to verify the well-known fact that a\ = —Ajv = —uj and that, 
if ai = a 2 = . . . = in Eq. ([]]), then /3i = j3 2 = • • • = 0. Eq. (Q) clearly shows the presence of 
several nonanalytic terms with exponents depending on 1/A, Aj/A, and 7/A. 

As pointed out by Alan Sokal, |2(], the nonanalyticity of the RG functions can also 
be understood within Wilson's RG approach. We repeat here his argument. Consider the 
Gaussian fixed point which, for 3 < d < 4, has a two-dimensional unstable manifold Ai u : The 
two unstable directions correspond to the interactions 2 and <p A . Then, notice that continuum 
field theories are in a one-to-one correspondence with Hamiltonians on M. u and that the FT 
RG is nothing but Wilson's RG restricted to M. u . But now Ai u has no special status at 
the nontrivial fixed point. In particular, there is no reason why it should approach it along 
the leading irrelevant direction. Barring miracles, the approach should be along a generic 
direction which has nonzero components along any of the irrelevant directions. If this happens, 
nonanalytic terms are present in any RG function. 

In order to clarify the issue, Ref. || determined the asymptotic behavior of (3(g) for g — > g* 
in the continuum field theory for N — > 00 and 2 < d < 4, and showed that Eq. (|2]) holds 
and that the expected nonanalytic terms are indeed present. In Ref. [24| the computation was 
extended to two dimensions, finding again nonanalytic terms. 

The presence of nonanalyticities gives rise to systematic deviations in FT estimates, as it 
did in MC and HT studies. In the next Section we will present a test-case and we will discuss 
the type of deviations one should expect. 
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2 A simple example 



In order to understand the role of the nonanalytic terms in the analyses of the FT perturbative 
expansions we have considered a simple zero dimensional example. Define 



/OO fjq- 
-= x* exp 
-oo V Z7T 



X 



gx_ 

24 



(3) 



For c and p not integer, this function has a branch point for g = 1 and thus it should mimic 
the behaviour we expect for FT expansions. For g — > 1, we have 



f(g; c,p) ^Z - Z x (l -g) + c(l - g) 1+p + 0((1 - gf)), 



(4) 



where Zq = 1.37556014, Z\ = —0.679325. We wish to repeat here the same steps performed in 
the calculation of g* and u. Therefore, we determine g* and Z x , by solving the equations: 

f(g*; c,p) = Z , Z x = /'(#*; c,p). (5) 

Of course, f(g;c,p) is replaced with an appropriate resummation of the its perturbative ex- 



pansion. We use here the resummation scheme proposed in Ref. |25| that makes explicit use 
of the location of the Borel-transform singularity, but similar results are obtained extending 
the Borel transform by means of Pade approximants (we mention that one could also use the 
perturbative series in the bare coupling [ESJ). The mean values and errors are determined by 
using the procedure of Ref. J27j. In the absence of nonanalytic term, i.e. for c = 0, using the 
nth-order expansion, we obtain 



n = 6: g* = 1.00025(131) Z x = -0.6791(178); 
n = 8: g* = 0.99997(10) Z x = -0.6800(18); 
n = 10: g* = 1.00000(1) Z x = -0.6791(2). 



There is good agreement, the precision increases by a factor of 10 every two orders, and the 
error bars are correct. 

Then, we consider the role of the nonanalytic corrections, by adding a term that is small 
compared to the analytic one. We choose c = —Z±/b. 
Now, for p = 1/10 we obtain 



n = 6: g* = 1.0043(62) Z x = -0.550(20); 
n = 8: g* = 1.0066(15) Z x = -0.550(3); 
n = 10: g* = 1.0062(5) Z x = -0.552(2). 



In this case the agreement is poor, especially for Z x , and, even worse, the errors are completely 
incorrect. This can be understood from Figs. |l| and ||] where we show the distribution of the 
approximants that are used. These distributions are nicely peaked, but unfortunately at an 
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Figure 2: Distribution of the results for the resummations of Z x for p = 1/10, c = —Z x /5. 



incorrect value of g* and Z\. Thus, in the presence of these (strong) nonanalyticities, the fact 
that the approximants have a narrow distribution is not a good indication that the results are 
reliable. Also, the stability of the results with the number of terms of the series is completely 
misleading. As we shall discuss below, this is what we believe is happening in two dimensions. 

If we instead consider a weak nonanalyticity, i.e. p ~ 1, the discrepancies we have found for 
p = 1/10 are much smaller, although still present. For instance, for p = 9/10 

n = 6: g* = 1.013(45) Z x = -0.755(270); 
n = 8: g* = 1.014(7) Z x = -0.655(23); 
n = 10: g* = 1.006(4) Z x = -0.668(2). 

In this case the results are consistent with the exact values, although the errors are still slightly 
underestimated. As expected, the largest discrepancies are observed for Z x . 

3 Conclusions 

In the previous Section we have shown that nonanalytic terms may give rise to systematic 
deviations and a systematic underestimate of the error bars. Now, what should we expect in 
the interesting two- and three-dimensional cases? 

In three dimensions A « 0.5 and A2/A is approximately |28| 2. Thus, the leading nonana- 
lytic term has exponent A2/A and is not very different from an analytic one. Thus, we expect 
small corrections and indeed the FT results are in substantial agreement with the estimates 
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obtained in MC and HT studies. However, small differences are observed for 7 and u for N = 
and N = 1. For instance (A = ujv) 



N = 7ft = 1.1596(20) Ref. 0, 7 M c = 1.1575(6) Ref. |29|, 

N = A FT = 0.478(10) Ref. 0, A MC = 0.517(7)1^° Ref. 

N — 1 7ft = 1.2396(13) Ref. §, 7ht = 1.2371(4) Ref. 0, 

j\T=l A FT = 0.504(8) Ref. 0, A MC = 0.533(6) Ref. BTJ. 



There are slight differences, especially for oj, but still at the level of a few error bars. Note that, 
as discussed in Ref. M, part of the error may be due to a slightly incorrect estimate of g*. 
Using the estimate of g* obtained from the analysis of the HT expansions, the FT estimates 
change towards the HT and MC values. 

Larger discrepancies are observed in two dimensions. For iV > 3 it is easy to predict 
the behavior of the RG functions at the critical point. Indeed, the theory is massive for all 
temperatures. The critical behavior is controlled by the zero-temperature Gaussian point and 
can be studied in perturbation theory in the corresponding iV-vector model. One finds only 
logarithmic corrections to the purely Gaussian behavior. It follows that the operators have 
dimensions that coincide with their naive (engineering) dimensions, apart from logarithmic 
multiplicative corrections related to the so-called anomalous dimensions. The leading irrelevant 
operator has dimension two EE 



33] and thus, for m — > 0, we expect 



g{m) 



g* |l + cm 2 In 



m 



l + O 



ln(— lnm 2 
lnm 2 



(6) 



where ( is an exponent related to the anomalous dimension of the leading irrelevant operator, 
and c is a constant. Therefore, 



(7) 



with Ag = g* — g. Clearly there are in this case logarithmic corrections and therefore, we expect 
large deviations in the determinations of to, which should be 2. These deviations are indeed 
observed: for iV = 3 the analysis of the five-loop series p|J yields the estimate (3 f (g*) = 1.33(2), 
which is very different from the expected result (3'{g*) = 2. 

For N = 1, we can repeat Nickel's analysis, using the fact that, from conformal field theory, 
we can compute the RG dimensions of all relevant and irrelevant operators. Using these results 
we predict thatQ 



>(i 



b^Ag^ 7 + b 2 \Ag\^' + hz \Agr< + 



2/7 



1 3/7 



3 It has been claimed sometimes that in the Ising universality class the leading irrelevant operator has u> = 4/3, 
so that P'(g*) — 4/3. This claim is incorrect. Indeed, such operator can only appear in nonunitary extensions 



of the Ising model, but not in the standard (unitary) 
Ref. \M. 



field theory. For a detailed discussion, see App. A of 
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where Ag = g* — g. Such an expansion is confirmed by an analysis of the lattice Ising model. 
Again we find strong nonanalyticities, and correspondingly we expect large deviations. And, 
indeed such large deviations are observed: The analysis of the five-loop series gives [[34] g* = 
15.39(25) and P'(g*) = 1.31(3), to be compared with the exact prediction /3'(g*) = 7/4 and the 
estimates g* = 14.69735(3) (Ref. JEJ) and g* = 14.6975(1) (Ref. |§). 
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